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We study the four-wave kinetic equation of weak turbulence linearized around the 
Rayleigh-Jeans spectrum when the collision integral is associated with short-range 
interactions between non-relativistic bosonic quasiparticles. The technique used for 
the analysis of the stability is based on the properties of the Mellin transform of the 
kernel in the integral equation. We find that any perturbation of the Rayleigh-Jeans 
distribution evolves towards low momentum scales in such a form that, when i — > cxd, 
all the particles occupy a sphere of radius arbitrary small. 

I. INTRODUCTION 

Since early work that has addressed the kinetics of Bose-Einstein condensation even before 
the experimental realization of BEC in weakly interacting atomic gases [l|,S,S,H,li, the 
description of the growth of the condensate has been the subject of considerable attention 0, 

as, 0,1 

llOj. At the first stage of evolution, a description based on the Boltzmann kinetic 



equation is adequate. In the homogeneous case, it has the form 

dtn{pi,t) = ^ I \M\'^{2Tc)^S{pi+p2-p3-P4)'2nS{uJi + uj2-uj3-^^4.) 



P2P3P4 

X [(1 + ni)(l + n2)n3n4- 721^2(1 +7^3) (l + n4)], (1) 



where we are using to denote momentum integration, 

d^p 



(2) 



Here n{p, t) is the distribution function, the average density of particles with momentum p 
at time t. The square-amplitude is taken as |A^p = (Svram"^)^, with a the scattering length 
parameterizing the Fermi pseudopotential V{ri—r2) = AiTam~^6{ri—r2) appropriate to low- 

2 

energy scattering by a interaction of finite range a. The dispersion law is Lu(p) = 

Based on the assumption that the formation of a condensate is a low energy process where 
the occupation numbers obey the condition n{t,p) ^ 1, it is customary to approximate the 



2 



Boltzmann equation by the kinetic equation of wave turbulence (ill . 



1 / Sna^ ^ 



dtui = - ( J / (27r) 5{pi +P2-P3- P4)27r5(u;i + iU2 - lus - u^) 

X [{rii + n2)n3n4 - nin2{ns + rii)] . (3) 



Using considerations of homogeneity [13[ , it is possible to obtain stationary power-law solu- 
tions n = Bp~'^, for the pair of values z/ = 3, 7/3. These Kolmogorov spectra are conjugate to 
the obvious solutions for z/ = 2, 0. However, the results of several studies which 
have addressed the question of the dynamics underlying Eq. seem to exclude any signifi- 
cant role of the Kolmogorov spectra in the formation of a finite time singularity correspond- 
ing to a condensate. Within the framework of kinetic theory, the dynamics before blow-up 
is rather described by a self-similar solution of the form n{t,p) = /^"^/^^((t* — 
where t^, is the time of blow-up and the parameters a and f3 are related by/5 = a — 1/2 (for 



a review, see Ref. [IJl)- The determination of the ratio a/ (3 poses a very difficult nonlinear 
eigenvalue problem which can not be solved by scaling arguments. The numerical observed 
value a//3 ~ 1.234 is significant different from the exponents z//2 = 3/2,7/6,1 associated 
with Kolmogorov and equilibrium solutions. However, the classical Ray leigh- Jeans spectrum 
with zero chemical potential plays a special role as boundary in the appropriate parameter 
space between the regions of condensate and normal phases (see Fig. [T]). To understand 
its significance it can be of some interest to study in detail the stability properties of this 
equilibrium solution. 

In this paper, we address the question of the linear stability of the Rayleigh- Jeans spec- 
trum {u = 2) of ([3]) whose perturbations could be related to processes driving particle 
transport towards the region of p = 0. Let us remark that if one tries to use the equilibrium 
n = T/uj{p) in the entire interval of time, one faces a serious difficulty. Consider a prescribed 
total particle density n to be accommodated by a singular distribution nQ(T){2n)^6{p) and 
an equilibrium Rayleigh- Jeans distribution through 

n = no{T)+ / =noT + — A. (4) 

Jq py{2m) (27r)^ 

For a given cut-off A, the second term grows with T, and when it reaches the prescribed 
total density n, the density of particles in the p = state vanishes. The temperature above 
which no = is Tc = 7r^n/(mA). In this model, the particle density in the condensate is 
given by 

no(T) =n(l-r/Te), T < T,. (5) 

From this fact, one would expect that the Rayleigh- Jeans spectrum can not be stable with 
respect to changes of temperature: a small change in T in the neighborhood of may 
move the system toward a state where no 7^ or not. Remarkably, this will be reflected 
by the fact that the linearized kinetic equation about n° = T/uj{jp) does not admit as 
a solution the function A{p,t) = 6T{p,t)/T = constant. This argument pointing out 
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an instability does not apply to the equilibrium distribution function, n oc p ^, which 
arises as the solution of the three-wave kinetic equation at the end of the condensation 



process la, ll7| • This comes from to the low frequency limit of the Bogoliubov dispersion 



law (jJsip) = \/4:Trap'^ /m + {p'^/{2m)y. In that case, the distribution function n corresponds 
to Goldstone-type excitations whose number is not conserved. Our focus will be on the 
regime in which the four- wave interactions are the most important. The transition from 
this stage to the three-wave regime has been recently studied within the two-dimensional 



Gross-Pitaevskii model in |18l. Il9l|. These authors have shown that this transition involves 
an intermediate step where strong interactions between topological defects take place in a 
similar way to the Kibble-Zurek mechanism in continuous phase transitions. 

It is important to point out that the evolution of an initial distribution function n{p, t = 0) 
in the regime of large occupation numbers produces a change in the entropy density 



S f 

— = / [(1 + n) ln(l + n) — nln 
^ J p 



n] (6) 



given by 

d{S/V) 1 fSira 



dt 4 V m 



(27r)^5(pi +P2-P3- P4)2vr5(a;i + 0)2-1^3- (^a) 
P1P2P3P4 



/ 1 1 1 M 

xnin2n3n4 \ , (7j 

V^i n2 ng Hi J 

which only vanishes for distribution functions with the shape of Rayleigh- Jeans. This result 
and the conservation laws imply that the final state of the evolution is consistent with the 
presence of a regular Rayleigh- Jeans distribution with or without a singular distribution, 
depending on whether the particle and energy densities in the initial distribution are below 
or above their critical values (see Fig. [T]). It follows that a perturbed Rayleigh- Jeans dis- 
tribution T/uj{p) + 6n{p,t = 0) can not evolve to the same = T/uj{p) unless the energy 
5{E/V) and the particle number 5{N/V) of Sn{p,t = 0) vanish. 

Although the theory of wave turbulence has been used in a variety of systems, ranging 
from ocean or capillary waves to plasma waves or elastic waves (see, for example, [20i] for 
a short review), specific studies on the stability of the stationary spectra are rather scarce. 
In this respect, another objective of this work is to present a specific new application where 



the techniques of Ref. [12j are used. 

This paper is organized as follows. In Sec. |IT1 we derive the linearized kinetic equation 
about the Rayleigh- Jeans spectrum. In Sec. IIIIl we introduce the Mellin function, and 
present some of its properties. The character of the evolution of the perturbation is con- 
sidered as p — > 0,00 and as t — 00. We also compute the leading asymptotics as t — 00 
of the perturbed particle number which is contained in the interval < p < A. Section IIVI 
gives our conclusions. Details of some calculations regarding the reduction of the integral 
equation to the Fredholm form and the evaluation of the Mellin function can be found in a 
pair of appendices. 
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FIG. 1: (color online). Parameter space at equilibrium for large occupation numbers. The curves 
are lines of /-i/T constant, obtained as parametric representations of the scaled quantities N/V 
and E/V in terms of the parameter A. The thick line corresponds to /i = 0. Each point on this 
curve corresponds to a Rayleigh-Jeans distribution with a critical temperature Tc = n'^N/ (mVA) = 
6'k'^E/{VA^) where N/V and E/V are the particle and energy density at the prescribed equilibrium. 
Each point in the region above the /i = line is associated with a pair (T, ^ < 0) of the equilibrium 
distribution n{p) = T/{uj{p) — ^u)). The dotted and dashed curves correspond to /i/T = —0.1 
and —0.4 respectively. For a point below the /i = line, the departure of its abcise from the 
curve corresponds to the particle density in the condensate and the energy density is given by its 
ordinate. 



II. THE LINEARIZED KINETIC EQUATION 



Here we derive the linearized kinetic equation about the Rayleigh-Jeans spectrum 

n\p) = Tuj{p)-\ (8) 



where it is understood that the temperature is at the critical value T 
in this model. We shall closely follow the conventions of 12 



= Tx^n/{mK) 
The dimensionless quantity 
A{p,t) represents the relative change in the particle number and may be either positive or 
negative. Substituting the expression 



n(p,t) =nO(p)(l + A(p,t)) 



(9) 
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into ([3]) and retaining the terms of 0{A) gives an integral equation for A{p,t), 

dtA = L[A], (10) 

with the result 

Z \ m J Jp2P3P4 

-n^{ps)n\p4)A{p2,t) + 2n\p2)n\p4)A{ps,t)} , (11) 
where we have used the conservation of energy 

1 1 _ 1 1 

n^{pi) rfi{p2) n^ips) n^ip^)' 

and the symmetry of the integrand under ps ^ p^ to duplicate the last term. As we 
shall see, one must be cautious using this rearrangement, which assumes regularity of the 
integrand. Here we require that the solution A{p,t) falls off sufficiently rapid as p ^ cxo 
to permit this kind of reordering in the integration. In the present case the several integral 
terms of the operator L containing the function A{p, t) for different arguments p are well 
defined separately, and it is possible to write 

L[A] = NA{p,t) + j U{p,p')A{p',t)d'p', (13) 

where the coefficient turns out to be a constant. 

Due to the conservation laws, one expects that the collision integral L[A] vanishes for 
certain forms of A. One of them, which reflects conservation of particle number, is A{p) oc 
This collision invariant corresponds to the fact that the kinetic equation ([3]) admits the 
Rayleigh- Jeans spectrum with non-zero chemical potential as an exact solution 

T 

<p) = , (14) 

uj{p) - /X 

which implies that A = fi/uj{p) is a zero mode of the linearized equation and 

^ = -/ ^Uip,p')d'p' (15) 

would be true. In fact, this turns out to be the case, as we have explicitly checked in 
Eq. flA24l) of the appendix. Since the linear momentum is conserved, a possible solution of 
Eq. da]) is 

T 

nip) = , (16) 

u{p) — V ■ p 

which would imply that A oc v ■ p is also a zero mode of the operator L. 
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Finally, from the conservation of the kinetic energy in collisions one would expect a zero 
mode A = constant but, surprisingly, this does not occur, at least for the operator L in the 
form ( ITTi) . Consider L[l] as follows from Eq. ( ITTi) : 

^[1] = ^ ( — — ) / {2irf6{pi+p2-p3-P4)'2n6{uji + U2-uj3-uJ4) 

X [-n°(p2)n°(p3) + n°(p2)n°(p4)] • (17) 

It it were allowable to make the exchange ps p^ in the second term, we would obtain the 
expected result L[l] = 0. However, one finds that 

/ 6{pi + P2 - P3 - P4)S{UJ1 +UJ2-i^3- C^4)?^°(P2)^°(P4) 
J P2P4 

'ln--l , P3-^oo, (18) 



271^ P3 V P3 



and 



/ S{pi + P2 - P3 - P4)S{UJI +UJ2-UJ3- C^4)^°(P2)n°(p3) 
Jp2P4 

'In 1 , p2 ^ 00. (19) 



27r^ pI \ p2 

These asymptotics give rise to leading divergences of the form ln^(A/pi) for the remanining 
integration over p^ and p2 respectively, being A an ultraviolet cutoff. This points out that 
it is not legitimate to use the previous rearrangement. We can see, however, that there is a 
cancellation between both divergences with the result that L[l] > 0. Using the notation in 
Appendix |A] we obtain 

/"°° dx Stt^ 
tL[1] =-fN+ Vo— = 7jv + 2G + — ^ 0.8225. (20) 
Jo X 8 

It is convenient to expand the kernel U{p, p') in terms of Legendre polynomials of the 
angle dpp> =p - p', 

00 

U{p,p') = Y,^i{p,p')Pi{cosepp>), (21) 

(=0 

and to expand [l^ the perturbation A in spherical harmonics relative to the orientation of 
P, 

A{p,t) = Y,AUP,t)Yr{p). (22) 

l,m 

From the addition theorem for spherical harmonics 

^/(cos V) = ^ E ^r(p)>^™(pO*, (23) 

m=—l 
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and the orthogonality relation {dQp' denotes the corresponding surface element of the three- 
dimensional unit sphere) 

j dnp,Yrip')Yr\pT = 5u'5mm', (24) 

a set of uncoupled equations emerges, each one of them labeled by (/, m): 

dtAira{p, t) = N{p)Aim{p, t) + J Ui{p,p')AUp\ t)p'^dp' . (25) 

In the homogeneous case, the angular momentum projection m may be suppresed because 
of isotropy, and it is enough to use / to label each irreducible perturbation. 

Since we are only interested in the / = mode, we drop this label to simplify the notation. 
In order to obtain the explicit form of the integral equation we must first evaluate the 
integrals involved in Eq. flTT]) . Some details of these calculations are given in Appendix 1X1 
One finds 

rd,A{p,t)=^^A{p,t)+ I^^V A{p',t)^, (26) 

where the dimensionless constant 7iv and r with dimension of time are given in (lASP and 
( 1A21I) . respectively. The kernel V(x) is given in Eqs. ( 1A20I) . ( lAllI) . ( 1A17I) and is depicted in 
Fig. [21 This equation is a special case of the general equation which arises in the linearization 
of the kinetic equation about the Kolmogorov spectrum. The resulting equation has the 
homogeneous form [l2| 

d^AM=p-^ j^U A{p',t)^. (27) 
In the present case, h = and U(x) is expressed as 7Ar5(a; — 1) + V(x). 



III. THE MELLIN FUNCTION AND THE ASYMPTOTIC BEHAVIOR OF 

THE SOLUTION 

The Mellin transform is a very useful tool for solving the integral equation fl26l) . In terms 
of the Mellin transform of V, we define 

POO 

W{s)=^jv+ Vix)x'-^dx = ^N + Wv{s), (28) 
Jo 

and thus the Mellin image J-'{s, t) of A{p, t) is a solution to 

TdtJ^{s,t) = W{s)T{s,t). (29) 

Hereafter we will use t and p to denote the dimensionless variables t/r and pi, where is 
some momentum scale entering into the initial condition A{p,t = 0). With these notation, 
the solution of the integral equation may be written as 

A{p,t) = — :/ p-'^{s)exp{W{s)t)ds, (30) 

271"^ J a—ioo 
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FIG. 2: (color online). A plot of the kernel V(x) for < x < 3. Near the singularity at x = 1, it 
behaves as 7r(2|x — 



where is the Mellin transform of the initial condition A{p,t = 0) 

POO 

= ^(s, t = 0)= / Aip, t = 0) p'-^dp, (31) 
Jo 

and a stands somewhere in —2 < a < 5, corresponding to the fundamental strip where W{s) 
is analytic. This property follows from the behavior of V{x) 

V(x) ~ -4x2 l^ln X - ^ + . . , X ^ 0, (32) 

4 / 2 \ 

V(x) ~ — (lnx + - + ...), X oo. (33) 



\ '3 

which shows that the Mellin function has double poles at s = — 2 and s = 5, 

"'t-') ~ ^ - (34) 

~ (7^. » - 5. (35) 

Some indications about the computation of the Mellin function are given in appendix 
iBl One finds that W{s) has zeros at s = 1 and s = 2 (see Fig. [3]), with the corresponding 
solutions A = p^^jP^"^. Another property is the symmetry with respect to the line Res = 
3/2, 

W{3/2 + s) = W{3/2- s). (36) 
The Taylor series about s = 3/2 reads 

W{s) 0.0786 + 0.3061(5-3/2)2, s 3/2. (37) 
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FIG. 3: (color online). A plot of the Mellin function W{s) for real values of s. The symmetry with 
respect to s = 3/2 is clear. The intercept is W{0) « 0.8225. 



A. The asymptotics as t | lnp| ^ +oo 



We begin the study of the stability by deriving from ( |30l) the infrared behavior of an 
evolved initial perturbation for arbitrary t. We require that a < 1 for the integral corre- 
sponding to the departure of the particle density 

Sn{t) = fn'{p)A{p,t)^^ (38) 

to converge as p 0. To evaluate (l30l) . we apply the steepest descent method to the integral 

A{p, t) = — I <l>{s) exp {xs + W{s)t) ds, (39) 

where x = —\np with x —>■ oo, and C is a contour which runs along a straight line parallel 
to the imaginary-s axis, with —2 < Res < 1. We assume that the initial perturbation has 
compact support, and its Mellin image vanishes sufficiently rapid as |Ims| oo to permit 
the existence of this integral. Note in passing that this property enables us to put the upper 
limit of the integral (!38l) as infinity. The saddle point occurs when W'{s) = —x/t. To find 
this location as x ^ oo we may use the asymptotic estimation and obtain 

s = -2 + 2 - , (40) 



which depends on x/t. A change of variable s = —2 + 2{t/xY^^u fixes the movable saddle 
point at u = 1, and the integral ( !39l) can be approximated by 



X 



1/3 



/ $(«) exp {x^/H'/%2u + 1/u^)) — , (41) 
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where = $(s('u)). Now the steepest descent method directly produces the asymptotics 
>4(p,t)~^^^^exp(3t^/=^(-lnpf/3), rMnp^-oo. (42) 

This resuh shows that, as t ^ +00, the initial perturbation becomes very large in the region 
where p > and t~^lnp — > —00, but preserves the infrared convergence of the integral for 
the particle density. 

A similar estimation using the steepest decent contour in the neighborhood of s = 5 — 
2(t/lnp)^/^ gives the leading behavior as Inp 00: 

AP,t) - ,.f\,,, exp (3t^/^(lnp)^/^), t-Mnp^oo. (43) 



where now a > 1 in Eq. fl30|) . for the integral 

(n f) - 

(27r) 

to converge as p — 00. 



B. The asymptotics as t ^| Inp \ — > 
In this regime, the relevant saddle point occurs near sq = 3/2 at 

Thus, the steepest descent evaluation yields the exponentially decreasing asymptotics 

A(p,t) ^iE2^= exp (W(3/2)t) , rHnp^O. (46) 

p3/2^27rW^"(3/2)t ^ ^ ^ 



The behavior of the function A{p, t) in this regime bears some resemblance with what is 



called interval instability in Notice nevertheless that the instability results in p^] can 



not be directly applied in our case since h = 0. 



C. The asymptotics of 5n{t;pQ) 

A more conclusive indication about the nature of the instability can be easily obtained 
by studying the induced change of the particle number density Sn{t;pQ) corresponding to 
the modes within the sphere p < Pq. Here Pq refers to an arbitrary scale which is smaller 
than the cut-off scale corresponding to the critical temperature. Let us consider (l38l) with 
a finite upper limit of integration p^. A first integration over p gives 

Snit;po) = ^ / ^<|.(.)e*^W^, -2 < Res < 1, (47) 
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where C runs parallel to the imaginary-s axis. Since the integrand falls off sufficiently rapidly 
at large |Ims|, the theorem of residues yields the exact expression 

6n{t-po) = ^$(l)e*^« + ^ / ^$(.)e*^(^)^, (48) 

TT^ TT^ J (J, 1 — S ITll 

where the new contour C runs parallel to the imaginary-s axis with 1 < Res < 5. The first 
piece is the contribution of the single pole at s = 1, which has been encircled in a negative, 
clockwise sense. Noting that W{1) = and making use of 

M0;oo) = — / n^{p)A{p,t = 0)p^dp= —^1), (49) 
In Jq vr 

the first piece of Eq. becomes time-independent and insensitive to the size posed by 
Pq. It corresponds exactly to the particle density coming from all the modes of the initial 
perturbation, and gives the leading behavior as t — oo. With regard to second term, the 
steepest descent evaluation along Res = 3/2 yields 

^$(s)e*^(^) "^(3/2) ^^(3/2), ^ ^ 

which shows an exponential decay law since W{3/2) ^ —0.0786 < 0. These results reveal 
that an arbitrary initial perturbation with support in some interval to the right of po, evolves 
towards the infrared momentum region. This gives rise to a situation where all the initial 
excess or defect of particles, depending on the sign of 6n{0,oo), appears finally located in 
the interval {0,Pq), irrespective of the magnitude of po- This is rather peculiar and seems to 
be consistent with some kind of generation of a large-scale structure. 



IV. CONCLUSION 

In this paper we have considered the problem of the stability of the Rayleigh- Jeans 
equilibrium solution of the theory of weak wave turbulence, in the case of the four-wave 
contact interactions with a non-relativistic dispersion law. It is believed that this theory 
captures the basics facts of the non-equililibrium processes which occur in weakly interacting 
Bose gases in the 'classical' regime where the occupation number is large. 

Based on the general theory of the stability of weak-turbulence Kolmogorov spectra 



presented in Ref. [l2j, we have first derived the explicit expression for the linear evolution 
equation for the scalar (/ = 0) perturbations, and then we have applied the method of Mellin 
function to analyze the stability of solutions. Had the index h been positive, the sufficient 
condition of Balk and Zakharov, W{0) > 0, would immediately imply the instability of the 



initial perturbation in the sense of Ref. [12]. However, this effective criterion for checking 
the stability is inapplicable in the present case. In order to see the the character of the 
evolution of an arbitrary initial perturbation, we have evaluated the leading asymptotic 
behavior as | Inp | oo with t held fixed, and as t ^ oo with p held fixed. But a better 
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understanding about the time evolution of a generic solution comes from the asymptotics for 
the integrated particle number as t — oo. This evaluation shows clearly the main feature of 
the evolution of the perturbation: the initial departure of the particle density concentrates as 
t — i> cx) into an interval of arbitrarily small width around the origin in momentum space. Let 
us remark again that numerical work based on the nonlinear kinetic equation strongly 
indicates a finite time for the formation of a singularity, according to the scenario described 
m Refs. 0, H, H, Esl where the nonlinear effects are important. We must then be very 
cautious in seeing the asymptotic instability of the hnearized equation around a Rayleigh- 
Jeans spectrum as a hint of the Bose-Einstein condensation. 

Finally, it is worthy mentioning that natural extensions of this work include the analysis 
of the stability of the Rayleigh- Jeans spectrum of three-wave kinetic equations such as those 
that arise at the final stage of the condensation process, and the stability for T > where 
the chemical potential is finite. Work in this direction is in progress. 
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APPENDIX A: REDUCTION TO FREDHOLM FORM 



Here we give some details about the reduction of the linearized kinetic equation for / = 
to Fredholm form. First we compute the constant N. The angular integration over the 
orientation of ps in the terms of the integrand of (ITT]) containing the functions A{pi,t) and 
A{p2, t) may be expressed in a series of Legendre polynomials of cos 6*12 = Pi ■ P2 



^ / z=o 



2™ efi_(p|Hj^, |,v;,p,,„,(,,), (Al) 



where denotes the Heaviside step function. Only the /o coefficient will be needed in order 
to compute and to derive the integral equation for ^0 By using the standard formula 

fi = ^ J Jimiu), (A2) 
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for the coefficients in the expansion f{u) = J2i fiPii.'^)i o^i^ finds for pi < p2 



27rm^, , 27rm ^ ^-, 

— e{pi - ps) H e(p3 - Pi) - e(p3 - ^2) 

P1P2 P2P3 



2'Km\/p\ + P2 ~ P3 
P1P2P3 



0(P3 - P2) - 0(P3 - \ pI + pI 



(A3) 



and for pi > p2 



(Pl>P2) 





(P3 



-Q[P2 - P3) H [0(P3 - P2) - 6(P3 - Pi)\ 



P1P2 



P1P3 



2'Km^p\ + P2 — pi 
P1P2P3 



0(P3 - Pi) - 0(P3 - \ pI + pI 



(A4) 



Note that foips) vanishes for p^ > \/ pf + pi- To derive the asymptotic relation (fTOll we may 
insert the expression of fo^ps) which produces 



n^{p2)n^{p?,)h{p2)pldp3 



pI 



ln^-1 

P2 



P2 00. 



(A5) 



With the notation p4 = \/pi+ pi. — pI, the integral over p^ affecting Ao^pi^t) adopts the 
form 



[-n^{p2)n^{P3) -n^{p2)n^{Pi) + n^{P'i)n^{Pi)\ fo{P3)pldp3 



P2 



.P2 



(A6) 



where 



Q{x) 



2x'^ In X 
l + x2 



2xlnx 



2x 



=arctanh 



x^ 



x 



VTT 



x^ 



2x 



:arctanh 



6(1 -x) 

e(x - 1). 



(A7) 



Thus the resulting three-dimensional p2 integral giving is independent of pi. The explicit 
computation shows that the radial integration of Eq. (1A6P is given by 



7n 



r Q(x)- = -2G - ^ - ^ In (3 - 2V2) In (3 + 2v^) 

•/ ^ 

^ In' (17 + I2V2) - Li2 (3 - 2^2) + 4Li2 (-1 + V2 
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-4.71057 



(A8) 



where G is the Catalan's constant and Li2 denotes the dilogarithmic function. Therefore, 
putting all factors together, we obtain the value of the constant A^ 



A = ma I , 



(A9) 



with dimension of 1/time. 
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A similar integration for the term containing AQ{p2,t) produces 



n {ps)n ip4)foiP3)P3dp3 



where 



R{x) 



2 In a; 



1 + X2 



X 



21nx 



x(l + x^ 



+ 



=arctanh 



x^ 



X 



x^ 



X 



=arctanh 



x^ 



,P2 



6(1 -x) 

Q{x - I) 



(AlO) 



x^ 



(All) 



The same procedure may be followed to evaluate the contribution of the coefficient of 
^o(P3;^) ill the integrand of flTTl) . Now the appropriate expansion in Legendre polynomials 

reads 



j 5{pi + P2-Pi- Pi)5{uJi + UJ2-UJ3- uji)£pid^p^ 



271171 



P2IP1 -PsI 



rO 1 



P2IP1 - Pal 



J](7,PKcos0i3), (A12) 



1=0 



where 



9'^r^'"\p2) 



+ P2 — P3 
P1P2P3 

271171 

H B(p2 -Ps), 



0(P2 - a/pI - ) - 0(P2 - Ps) 



P2P3 

2nm^, , 27rm^, , 
B(P3 -P2) H B(p2 -Paj- 



The asymptotic behavior of (70(^2) is 

27Tm 



27rm 



5-0(^2) 0(P3 - Pi) H 0(Pi - Ps), P2 ^ 00, 

P2P3 P1P2 



(A13) 
(A14) 

(A15) 



but it is not dangerous for the remaining p2 integration of n^{p2)n^{PA), which is seen to be 

Snin^T^ ^ fpi 



11 {p2)n {pA)go{p2)P2dp2 



where 



S{x) 



2 In a; 



+ 



xVl — X' 



+ 



1 — 
2 Inx 

X(x2 - 1) ' xy/. 



: arctan 



pI 



X 



(A16) 



x^ 



arctan 



x^ 



6(1 -x) 

0(x - 11 



x^ 



The behavior 



5(x) ^ -21nx + 2, 



(A17) 



(A18) 
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determines the asymptotics which has been used in Eq. (fT8i) . 

Assembhng the various contributions produces the integral equation 

TdtAip, t) = 7^^(p, t) + Tv A{p\ t)^, 

Jo \P J V 

where V(x) is given by 



TT 



32ma2T2' 

Near x = 1 the kernel exhibits the behavior 



(A19) 



V{x)=R{x)\S{x), (A20) 

and the constant r is 



(A21) 



V(x) ~ -j^=. ^ - 1- (A22) 

The kernel has the properties 

/ V(x)— = 2G+^^ 5.5330, (A23) 
Jo ^ ° 

poo 

/ V{x)xdx = -jN, (A24) 
Jo 

the second assuring that p^^ is a zero mode of L. This completes the derivation of 

the explicit expression for the integral equation. 



APPENDIX B: COMPUTATION OF THE MELLIN FUNCTION 

The integrals involved in the computation of the Mellin function containing logarithmic 
terms can be expressed in closed form in terms of the polygamma function ip^'^\z) for n = 1: 

^x'^lnx , 1 ,,/3 + a\ 1 + 

^^^,.^-l^'(l±^y Re„>-1. (B2) 

To perform the remainder integration involving the inverse trigonometric functions, it is 
convenient to use the following series expansion 

u ^ (2n)! 2 

arctanw = , > ^ , r— , u < oo, (-d3) 

yrr^ii^ 22"(n!)2(2n + l)(l+M2)n' v ; 



and its counterpart for arctanhu. The integration term by term of the resulting series and 
the subsequent summation produces some combination of hypergeometric functions 3F2 and 
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2-Fi evaluated at ±1. This procedure yields 



X 



; arctan 



X 



dx 



^ r (g/2) 

2 r(l/2 + s/2) 

1 1 s_ 3 1 + s 
2'2'2'2' 2 



xVl + 



arctanh 



X 



X' 



dx 



X3F2 
1 



(B4) 



s{i-s: 



„ , s s 1 + s s 
3F2 1,-,-;^— + 



2' 2' 



V21n(l + V2) 
s- 1 



2Fi(l,£;l±^;-l].(B5) 



The other integrals in the interval [1, 00] are obtained from these by making the change 
of variable x — > 1/x. Assembling the various pieces, one finds the integral for the Mellin 
transform of V{x) 



Wy{s) 



V{x)x'-^dx 

2v/21og {I + V2) 

s-1 
2V21og(l + V2) 



s-2 



.F,(l,^;i^;-l 



2 ' 2 ' 



r(^) 
^r(^) 



r(2-|) 

2 



3F2 



1 1 s 3 1+s 



2'2'2'2' 2 



113-53^ s.-^ 
2' 2' 2 '2' ~ 2' 



(s-l)s 

2 



s-3)(s-2) 
3-s 



^^^l^'2'2'^'^+2'~^ 



3F.(1,^,^;2-J,^;-1 



2 ' 2 
3 -s 



4^ V 2 



2y 



(B6) 



which converges in the strip — 2 < Res < 5. 
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